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Abstract 

We analyze the chiral restoration and deconfinement transitions in the framework of a non-local 
chiral quark model which includes terms leading to the quark wave function renormalization, and 
takes care of the effect of gauge interactions by coupling the quarks with the Polyakov loop. Non- 
local interactions are described by considering both a set of exponential form factors, and a set 
of form factors obtained from a fit to the mass and renormalization functions obtained in lattice 
calculations. 
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I. INTRODUCTION 

The detailed understanding of the behavior of strongly interacting matter under extreme 
conditions of temperature and/or density has become an issue of great interest in recent 
years. Unfortunately, even if a significant progress has been made on the development of ab 
initio calculations such as lattice QCD, these are not yet able to provide a full understand- 
ing of the QCD phase diagram and the related hadron properties, due to the well-known 
difficulties of dealing with small current quark masses and finite chemical potentials. Thus 
it is important to develop effective models that show consistency with lattice results and 
can be extrapolated into regions not accessible by lattice calculation techniques. Recently, 
models in which quark fields interact via local four point vertices and where the Polyakov 
loop is introduced to account for the confinement-deconfinement phase transition (so-called 
Po,ya.ov-Na.bu-WLa.Mo (PNJL) .ode. Q B. 3 fl 3) have .ece.ved eo^.de.able 
attention. Here, we consider a non-local extension of these PNJL models, which includes 
terms leading to the quark wave function renormalization. Two different parameterizations 
are used: an exponential form, and a parametrization based on a fit to the mass and renor- 
malization function obtained in lattice calculations. In the context of this type of model the 
properties of the vacuum and meson sectors at T = /i = have been studied in Ref. 6|]. 

This contribution is organized as follows. In Sec. 2 we introduce the model lagrangian 
and its parameterizations. In Sec. 3 we present and discuss our results for the behavior of 
some thermodynamical properties and the corresponding phase diagrams. Finally, in Sec. 4 
our main conclusions are summarized. 



II. THE MODEL AND ITS PARAMETRIZATIONS 

We consider here a nonlocal SU(2) chiral quark model which includes quark couplings to 
the color gauge fields. The corresponding Euclidean effective action is given by 



G 

d'^x ^'ip{x)[-ip + m]ip{x)-^ jaix)jaix)-jp{x)jp{x) +W($[v4(x)]) 



(1) 



where ip is the Nf = 2 fermion doublet ip = {u,d)'^, and rh = diag{mu,md) is the current 
quark mass matrix, in what follows we consider isospin symmetry, that is = m„ = m^. 



The fermion kinetic term includes a covariant derivative D^^ = d^ — iA^^ where are color 
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gauge fields. The nonlocal currents ja{x), jp{x) are given by 

Ja{x) = j d!^z g{z) ^{x + ^ Ta'ipix-^ , 

Jp(^) = / d'z f{z) + I) ^ - I) > (2) 

Here, Ta = (l,z75r) and u{x') d v{x) = u{x')dxv{x) — dx'u{x')v{x). The functions g{z) 
and f{z) in Eq.([2]), are nonlocal covariant form factors characterizing the corresponding in- 
teractions. The scalar-isoscalar component of the ja{x) current will generate the momentum 
dependent quark mass in the quark propagator, while the "momentum" current, jp{x), will 
be responsible for a momentum dependent wave function renormalization of this propagator. 

To proceed it is convenient to perform a standard bosonization of the theory. Thus, we 
introduce the bosonic fields <7i^2{x) and vra(x), and integrate out the quark fields. In what 
follows, we work within the mean- field approximation (MFA), in which these bosonic fields 
are replaced by their vacuum expectation values o"i^2 and 7ta = 0. Next, we extend the 
so obtained bosonized effective MFA action to finite temperature T and chemical potential 
fi using the Matsubara formalism. Concerning the gluon fields we will assume that they 
provide a constant background color field ^4 = iAq = ig 6f^oG^X"'/2, where are the 
SU(3) color gauge fields. Then the traced Polyakov loop, which is taken as order parameter 
of confinement, is given by $ = iTr exp(i/50), where f3 = 1/T, = iAo. We will work 
in the so-called Polyakov gauge, in which the matrix is given a diagonal representation 
(j) = ^sAs+^gAg, which leaves only two independent variables, 03 and 08- Owing to the charge 
conjugation properties of the QCD Lagrangian [7], the mean field value of the Polyakov loop 
field $ is expected to be a real quantity. In addition, we assume as usual that 03 and 08 are 



real-valued fields [s], this implies that 08 = 0, then $ = [2cos(03/T) + l]/3. 

Within this framework the mean field thermodynamical potential results 



^2 rr2 



+ ^+m.T), (3) 
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Here, the shorthand notation Jp^ = J2n I d^vK'^'^Y has been used, and M{p) and Z{p) are 
given by 

M{p) = Z{p) [rrif + gip)] , Zip) = [1 - f{p)]-' , (4) 
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where g{p) and f{p) are the Fourier transform of g{z) and f{z), respectively. In addition, 
we have defined 



{2n + l)nT -ifx + 



(5) 



where the quantities 0c are given by the relation = diag{(j)r, 4>g, 4>b)- Namely, 0c = c 03 
with c = 1,-1,0 for r,g,b respectively. At this stage we need to specify the explicit form of 
the Polyakov loop effective potential. Here, we used the fit to QCD lattice results proposed 



in Ref . 
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fi'^^'^ turns out to be divergent and, thus, needs to be regularized. For this purpose we 
use the same prescription as in Ref. By- Namely 



OMFA 



MFA 



(6) 



where ^l^"^^^ is obtained from the first term in Eq. ([3]) by setting o"i = (T2 = and ^freg) 
regularized expression for the quark thermodynamical potential in the absence of fermion 
interactions. 



Qfree 
(reg) 



-4 T 



{2n) 



3 

c,k 



In 



l + e 



(7) 



with k = ±1. Finally, note that in Eq.([6]) we have included a constant Qq which is fixed by 
the condition that vanishes at T = /i = 0. 

The mean field values o"i 2 and 03 at a given temperature or chemical potential, are 
obtained from a set of three coupled "gap" equations. This set of equations follows from 
the minimization of the regularized thermodynamical potential, that is 



'MFA 



i-'i 'MFA 



^^MFA 







(8) 



d(7i da2 (903 

Once the mean field values are obtained, the (T, fi) behavior of other relevant quantities can 
be determined. 

In order to fully specify the model under consideration we have to fix the model parame- 
ters as well as the form factors g{q) and f{q) which characterize the non-local interactions. 
Following Ref. jo], we consider two different type of functional dependencies for these form 
factors. The first one corresponds to the often used exponential forms. 



f{q) = exp[-gVA?] . 



(9) 



4 



Note that the range (in momentum space) of the nonlocahty in each channel is determined 
by the parameters Aq and Ai, respectively. Fixing the T = fi = values of rric and chiral 
quark condensate to reasonable values = 5.7 MeV and {qq)^^^ = 240 MeV the rest of 
the parameters are determined so as to reproduce the empirical values /^r = 92.4 MeV and 
= 139 MeV, and Z{0) = 0.7 which is within the range of values suggested by recent 
lattice calculations . In what follows this choice of model parameters and form factors 
wiU be referred as parametrization SI. The second type of form factor functional forms we 
consider is given by 

9W = -r, 7TT ' W = ^ TT^fziq) , 10) 



where 



fM = 1 + {qVKY^" ' ; m = [1 + (gVA?)] ■ (11) 



5/2 



As shown in Ref.p'], with a convenient choice of parameters one can very well reproduce 
the momentum dependence of mass and the renormalization function obtained in a Landau 
gauge lattice calculation as well as the physical values of m^^ and In what follows 
this parametrization will be referred as S2. Finally, in order to compare with previous 
studies where the wavefunction renormalization of the quark propagator has been ignored 
we consider a third parametrization (S3). In such case we take Z{p) = 1 (setting f{p) = 0) 
and exponential parametrization for g{p). The values of the model parameters for each of 
the chosen parameterizations are summarized in Table I. 

TABLE I: Sets of parameters. 



SI S2 S3 



rric 


MeV 


5.70 


2.37 


5.78 


GsAl 




32.030 


20.818 


20.650 


Ao 


MeV 


814.42 


850.00 


752.20 


Kp 


GeV 


4.180 


6.034 




Ai 


MeV 


1034.5 


1400.0 
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III. RESULTS 



We start by analyzing the behavior of some mean field quantities as functions of T and 
Since the results obtained for our three different parameterizations are qualitatively quite 
similar we only present explicitly those corresponding to the parametrization SI. They are 
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FIG. 1: Mean fields ai, a2 and $ as functions of T for low (left), high (right) and CEP (central) 
chemical potentials. Note that the scale to the left corresponds to that of ai while that to the right 
to a2 and Since a2 turns out to be negative we plot — a2. 

given in Fig.l where we plot ai , (J2 and <!> as functions of T for some values of the chemical 
potential. Fig. la shows that at = there is a certain value of T at which ai drops rapidly 
signalling the existence of a chiral symmetry restoration crossover transition, its position 
being determined by the peak of the chiral susceptibility. At basically the same temperature 
the Polyakov loop l» increases which can be interpreted as the onset of the deconfinement 
transition. As /i increases there is a certain value of /i = ficEP above which the transition 
starts to be discontinuous. At this precise chemical potential the transition is of second 
order. This situation is illustrated in Fig.lb. The corresponding values {Tcep, /J'Cep) define 
the position of the so-called "critical end point". As displays in Fig.lc, for // > /jLcep the 
transition becomes discontinuous, i.e. of first order. Finally, for chemical potentials above 
IJ,c{T = 0) ~ 310 MeV the system is in the chirally restored phase for all values of the 
temperature. It is important to note that although cx2 appears to be rather constant in 
Fig.l, at higher values of T it does go to zero as expected. Concerning the deconfinement 
transition we see that as /i increases there appears a region where system remains in its 
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symmetry has been 
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confined phase (signalled by $ smaller than ~ 0.3) even though chira^ 
restored. This corresponds to the recently proposed quarkyonic phase 

The phase diagrams corresponding to our three different parameterizations are given in 
Figl2i Here the dotted line corresponds to the line of crossover chiral transition while the 
full line to the line of first order chiral transition. The dashed lines correspond to the 
deconfinement transition (the lower and upper lines correspond to $ = 0.3 and $ = 0.5, 
respectively). Comparing those of SI and S3 we see that the main effect of the wave function 
renormalization term is to shift the location of the CEP towards lower values of T and higher 
values /i. Concerning the lattice adjusted parametrization S2 we observe that it leads to 
even lower values of Tqep and higher values [Icep- 

250 

200 




100 200 



100 200 300 
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100 200 



FIG. 2: Phase diagrams for the three parameterizations considered. SI and S2 include quark wave 
function renormalization while S3 does not. SI and S3 correspond to exponential form factors while 
S2 to lattice motivated form factors. The dotted line corresponds to the line of crossover chiral 
transition and the full line to that of first order chiral transition. The dashed lines correspond to the 
deconfinement transition (the lower and upper lines being for $ = 0.3 and ^ = 0.5, respectively). 



IV. SUMMARY AND CONCLUSIONS 

A non-local extension of the PNJL model momentum which leads to momentum depen- 
dent quark mass and wave function renormalization has been studied. This model provides 
a simultaneous description for the deconfinement and chiral phase transition. Non-local 
interactions have been described by considering both a set of exponential form factors, and 
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a set of form factors obtained from a fit to the mass and renormalization functions obtained 
in lattice calculations. The resulting phase diagrams turn out to be qualitative similar, the 
position of the critical end point being the feature which depends more crucially on each 
particular parametrization. 
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